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ABSTRACT 

Over the last 15 years, around a hundred very young stars have been observed in the central parsec 
of our Galaxy. While the presence of young stars forming one or two stellar disks at « 0.1 pc from the 
supermassive black hole (SMBH) can be understood through star formation in accretion disks, the 
origin of the S stars observed a factor of 10 closer to the SMBH has remained a major puzzle. Here 
we show the S stars to be a natural consequence of dynamical interaction of two stellar disks at larger 
radii. Due to precession and Kozai interaction, individual stars achieve extremely high eccentricities 
at random orientation. Stellar binaries on such eccentric orbits are disrupted due to close passages 
near the SMBH, leaving behind a single S star on a much tighter orbit. The remaining star may be 
ejected from the vicinity of the SMBH, thus simultaneously providing an explanation for the observed 
hypervelocity stars in the Milky Way halo. 

Subject headings: black hole physics — Galaxy: center — methods: n-body simulations — stellar 
dynamics 



1. INTRODUCTION 

The central parsec of our Galaxy harbors different 
groups of young stars orbiting the central SMBH. These 
include one or two stellar disks at k, 0.1 pc from Sg r 
A* made up of 5-6 Myr old stars (jPaumard et al.ll2006l ). 
and a number of stars on highly eccentric orbits even 
closer, the so-called S stars (Eckart & Genzcl 1997]). 
With orbital periods as small as 15 years, the S stars 
provide powerful constr aints on the mass and size of 
the central dark object (|Ghez et al.ll2005D . While the 
origin of the stellar disk(s ) ob served in the Galactic 
center can be und erstood (iLev in & Bcloborodov 200 
Genzel et all 120031: 'Navakshin fc Cuadraj2005|; iGerhari 



20011: iMcMillan fc Portegies Zwart .2005t) . the S stars 
show a puzzling comb ination of interesting features 
([Eckart k, Genzel Il997l) : They are very young (< 
10 Myr), with less than 0.01 pc a factor of 10 closer to 
the SMBH than the closest young population known (the 
stellar disks) , and move around Sgr A* on randomly ori- 
ented eccentric orbits. 

A number of mechanisms to create the S stars have 
been proposed, but so far none of them was able to 
account for all their properties (see also review by 
lAlexandei1l2OO50 : It has been suggested that the S stars 
are not young but rejuvenated stars or exotic objects, 
but this is in conflict with their spectral properties, 
showing they are ordinary main-sequence B-type stars 
(iGhez et al.l 12005": 'Eisenhau er et al.ll2005t iMartins etlHI 
|2008|). In situ formation is excluded by the inexplica- 
bly high pressure required to overcome the strong tidal 
force s, and will also lead to a stellar disk (jGhez et al.l 
|2005| ). which cannot account for the random orien- 
tation and hig h eccentricity of the orbit s within the 
stars' Hfetime ()Giirkan fc HopmanI |2007[ ). The lat- 
ter argument also argues against the sugge sted sce- 
nario of an infalling cluster of young stars (jGerhardI 
l2001f). which would moreove r require a very high mass 
( Portegies Zwart et al.l |2003[ ). or an intermediate-mass 
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black ho le in its center for it to survive the tidal force of 
Sgr A* (jHansen fc Milosavlie"\^l2003l). Exchange cap- 
ture with compact remnants ( Alexander fc Liviol l2004( l 
only works with a large cluster of black holes and fails 
to explain the youth and high initial eccentricities. The 
Lense-Thirring effect rnay ac count for the innermost star 
(jLevin fc Beloborodovll200^ . S2, being dragged out of 
the disk plane, but cannot explain the orbital orienta- 
tion of the longer-period S stars, nor the proximity of 
the orbit itself. Interaction of stars in a central cluster 
with a fossil accret ion disk can lead to h igh eccentricities 
and orbital decay (jSubr fc Karasl[2005l) if the disk mass 
is sufficiently high, but fails to explain the observed bias 
toward young B-type stars. Furthermore, such an ac- 
cretion disk is expected to fragment and form stars and 
should thus have been observed. 

Tidal disruption of binaries falling into the Galactic 
center on highly eccentric orbits has been shown to be 
an effective mechanism to create S stars, but until now 
requ ired very massive star clusters passing close to Sgr 
A* (iGould fc QuillenI [2003f) or the presence of a large 
enough population of compact massive perturbers in 
combination with many young stars on eccentric orbits 
(|Peretset al.ll2007[ ). 

In the following we show how the dynamical interaction 
of two stellar disks at larger radii, as they are observed 
in the Galactic center, naturally leads to the creation of 
S stars. Section [5] describes the setup of our numerical 
simulations. In § [3l we present the outcome of these 
simulations, showing how the disk stars can eventually 
gain very high eccentricities. Section |4] explains how bi- 
nary stars on such eccentric orbits around the SMBH get 
disrupted, leaving behind an S star and possibly a hyper- 
velocity star. In §[5] we discuss our results, showing how 
the mechanism suggested here explains all the observed 
properties of the S stars. 

2. MODELING DISK INTERACTION 

An extensive analysis of young massive stars in the in- 
nermost parsec of our Galaxy shows that almost all stars 
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Fig. 1. — Aitoff projection of normal vectors of the orbits of stars 
situated in two interacting disks. The separation of two points 
in the projected plot indicates the relative inclination of the two 
corresponding stellar orbits, and thus the disk thickness. The red 
squares indicate the initial state of two nearly flat disks inclined 
at an angle of 130° . The orientation of the disks changes with 
time due to precession. By 5Myr (circles), the massive disk shows 
a slight warping (bottom). The less massive disk (top to left) 
precesses faster and is affected more strongly by warping, especially 
towards the inner edge (green and yellow circles) , where it loses its 
disk shape. 

reside in two 6 ± 2 Myr old rotatin g disks — except for th e 
S stars in the central arcsecond (jPaumard et alJ [2006f ) . 
While the clockwise disk is very distinct, there is some 
debate over the disk s tructure of the counterclockwise 
system ()Lu et al.ll2006l ). However. [Cuadra et all (|2008D 
have shown that a single stellar disk cannot explain the 
high inclinations and eccentricities observed. Using di- 
rect iV-body calculations as described below, we find that 
the distribution of eccentricities and inclinations of the 
observed young stars' orbits is consistent with an initial 
configuration of two relatively flat disks, since the less 
massive disk is destroyed over time due to dynamical in- 
teraction (see §[3]). 

To investigate the stellar dynamical evolution of two 
interacting stellar disks, we have computed a num- 
ber of A'^-body integrations. The parameters of our 
model are chos e n so as to best fit the observations by 
iPaumard et al.l ()2006[ ) and are as follows: Our model 
consists of a 3.5 x 10^ Mq SMBH and two flat circular 
stellar disks, mutually inclined at an angle of 130°, with a 
surface density proflle that scales with distance as r~^-^. 
The disks have well-defined initial extents of 0.05 — 0.5 pc 
and 0.07 — 0.5 pc, respectively. 

We ass ume that the lower m ass limit of the stars ob- 
served bv IPaumard et al.l (|2006D is 20 M©. For our mod- 
els, we assume a disk age of 5 Myr, and hence an upper 
initial mass limit o f 42.5 M© for the o bserved disk stars 
(jHurlev et al.ll2000tlFuchs et al.ll2006f) . Using the initial 
mass function dN/dm ex m ~^'^^ as derived from the K- 
band luminosity function bv IPaumard et al.l l|200Q) with 
mass limits 1 and 120 M©, this leads to a total initial 
stellar mass of 10, 380 Mq to explain the 73 observed O 
stars. Allowing for 30% binaries and another 10% for 
stars disrupted by the SMBH or unobserved, and using 
a mass ratio of 2 : 1 for the disks, we have initial stel- 
lar masses of 9,900 and 4,95OM0, respectively. As a 
consequence, our disks initially consist of 535 and 267 
multimass stars (and binaries), respectively. The above 
parameters have proven to result in a distribution of ec- 
centricities and semimajor axes as well as an inclination 
between the two disks after 5 Myr which best match the 
observations. 



We used our bhint code (jLockmann fc Baumgardtl 
|2008[ ) for the integration, which has been developed 
specifically to calculate the dynamics of stars orbiting 
a SMBH, and includes post-Newtonian treatment up to 
order 2.5 to account for the effects of general relativ- 
ity. Furthermore, we have added the SSE package by 
iHurlev et al.l (|2000f) to account for the effects of stellar 
mass loss. 

The stars and binaries are modeled as point masses, 
and all binaries are assumed to be equal-mass. Tidal 
disruption is not considered; however, we assume that 
stars on very tight orbits (a < 80 AU) are swallowed by 
the SMBH. We switch on the post-Newtonian terms for 
the cen tral motion around the SMBH i n our bhint inte- 
grator (jLockmann fc Baumgardt) *2008l for eccentricities 
e > 0.9, and switch it off when e falls below 0.8. A fully 
post-Newtonian calculation shows that in the regime con- 
sidered in this Letter, relativistic effects are negligible 
for orbits with lower eccentricities within the considered 
timescale. 

We finally note that the detailed choice of parameters 
seems not to be too important. A number of models with 
variations in total mass, mass function, and disk age lead 
to comparable results. 

3. ACHIEVING VERY HIGH ECCENTRICITIES 

Due to mutual torques, the two stellar disks described 
above precess about each other, with a frequency propor- 
tional to the other disk's mass: The precession frequency 
of a star orbiting a SMBH of mass Msmbh at radius R at 
an inclination /3 relative to a narrow dis k of mass Mdi.sk at 
a radius i?disk can be approximated as (|Navakshinll2005f) 
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As can be seen, the precession frequency depends 
also on the dista nce of a star from the central SMBH 
()Navakshinll2005f ) . and since both disks have a finite ex- 
tent, stars at different central distances precess with dif- 
ferent frequencies, thus warping the originally fiat disks. 
In particular, the innermost stars quickly get inclined 
with respect to their disk of origin, and so start precess- 
ing about both disks. This leads to further thickening 
and deformation of the disks as well as large inclinations. 

Figure[l]depicts the evolution of the stellar disks' shape 
and thickness. Being exposed to the larger torque, the 
less massive disk is warped much faster, explaining why 
the observed thickness of the less massive disk in the 
center of our Galaxy is substa ntially higher, making its 
present existence controversial (jLu et al. 2006). 

Apart from pre cessioii, all stars are subject to the 
Kozai mechanism ()Kozailll962D : Stars on highly inclined 
orbits relative to an axisymmetric perturbation periodi- 
cally obtain small inclinations, and simultaneously gain 
high eccentricities. At their pericenter, stars on orbits 
with eccentricity e > 0.99 pass the central SMBH within 
~ 100 AU, less than the S stars' pericenter distances. 
The effects of general relativity causes eccentric orbits to 
precess, which may damp the Kozai mec hanism an d pre- 
vent further growth in eccentricity (Holm an et al.l [l997). 
However, our simulations show that in the presence of 
very massive, nearby disks, relativistic effects do not pre- 
vent eccentricities as high as e = 0.999. 
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Fig. 2. — Eccentricity evolution of stellar orbits. The top panel 
shows as an example the time evolution of orbital eccentricity of 
three stars with short Kozai periods. The cyclic nature of the 
Kozai mechanism, leading to extremely high eccentricities, is evi- 
dent. The solid and dotted lines in the bottom panel depict the 
fraction of stars having reached eccentricities higher than 0.9 and 
0.99, respectively, as a function of time. Any star on a Kozai oscil- 
lation keeps a high eccentricity for a short time only, but as most 
stars have a Kozai period longer than the computed time, these 
numbers are increasing until the end of the integration. 

Figure m shows the eccentricity evolution of some stars 
in our simulation, and the total fraction of stars that 
have reached eccentric orbits within a given time. A sig- 
nificant fraction of the stars obtain very high eccentric- 
ities. In contrast, a single stellar disk around a SMBH 
evolves only very slowly: In a corresponding integration, 
no star obtained an orbital eccentricity above 0.8 within 
a comparable time span. This is mainly due to the lack 
of an external torque, and the lack of large inclinations 
required for the Kozai mechanism to become effective. 

4. DISRUPTION OF STELLAR BINARIES 



As most stars form in binaries (jGoodwin et al.ll2007f ) , it 
is natural to assume that a significant fraction of stars in 
the stellar disks are followed by a companion. The orbits 
of these binaries around the SMBH are subject to the 
same Kozai mechanism as single stars, but a very close 
pericenter passage may break up the binary and leave 
one star on a tightly bound orbit around the SMBH. 

In a set of three-body integrations of a stellar binary 
orbiting a SMBH, we have investigated the distribution 
of binary disruption remnants. As can be seen in Figure 
[21 the eccentricity evolution due to the Kozai effect is 
smooth. 

We calculated a set of models including a SMBH, an 
intermediate-mass black hole (IMBH) to drive the Kozai, 
and a stellar binary on an initially circular orbit at high 
inclination with respect to the IMBH's orbit in order to 
test the impact of a Kozai oscillation on the binary's 
inner orbital parameters. Figure [3] shows the time evolu- 
tion of a 0.1 AU binary's both outer and inner eccentric- 
ity and semimajor axis. It can be seen that until binary 
disruption, only the outer eccentricity is changing. In 
particular, the binary does not widen. Only during the 
last few pericenter passages does the inner eccentricity 
grow to a value of e = 0.6, but this does not significantly 
change the resulting semimajor axis of the bound mem- 
ber, nor the velocity of the ejected star. Computations 
with wider binaries (1 AU) yield similar results. 

We thus conclude that any tight binary undergoing 
Kozai resonance reaches its point of disruption basically 
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Fig. 3. — Evolution of orbital parameters of a stellar binary 
undergoing Kozai resonance. The long-dashed and dotted lines 
show the outer and inner eccentricity of the binary, respectively, 
while the solid and short-dashed lines show the outer and inner 
semimajor axis, respectively. It is seen that the binary's inner 
orbital parameters stay constant until binary disruption, while the 
outer eccentricity grows due to the Kozai mechanism. 
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Fig. 4. — Orbital parameters about the SMBH of the (more 
strongly) bound member star of a disrupted equal-mass binary as 
a function of initial binary separation. While the disruption of wide 
binaries leaves the individual stars on orbits similar to the initial 
trajectory, the bound star of a disrupted tight binary ends up on a 
very tight orbit around the SMBH. The solid line depicts the initial 
semimajor axis before disruption. The shaded area indicates the 
semimajor axes of the S stars in the Galactic center. 

unperturbed. We have therefore computed a set of mod- 
els where the pericenter distance of a binary orbiting a 
SMBH is equal to the disruption distance (Fig. [4]). We 
find that the disruption remnants of binaries with sep- 
arations larger than a few AU orbit the central SMBH 
on similar trajectories as the initial binary, and so can 
still be observed as disk stars. However, the disruption 
of a tight binary (with separation of 1 AU or less) forces 
one star onto a very tight and eccentric S-star-type orbit, 
while the other star gains a significant amount of energy 
and may even be ejected from the Galaxy, thus becom- 
ing a progenitor of t he hypervelocity stars observed in the 
Gala ctic halo (e.g.. iBrown et al.]|2007l : [Gualandris et al.l 
|2005| ). For example, disruption of binaries with a separa- 
tion of 0.2 AU leaves one star on an orbit with semimajor 
axis consistent with those of the observed S stars, while 
the other star is ejected from the Galaxy with a velocity 
of up to 1000 km s~^ in more than half of our integra- 
tions. 

Binary disruption has been shown to be an effec- 
tive m echanism to create the S stars (|Gould fc Quilleia 
|2003[ ). One way to create young close-passage binaries 
is by massive perturbers like star clusters or molecular 
clouds at distances of a few pc from the Galactic center 



(jPerets et al.ll2007l) . Our simulations show that these bi- 
naries are also a natural consequence of the interaction of 
the two stellar disks: Using a binary fraction of 30%, 18 
of the B-type binaries achieve pericenter distances be- 
low their tidal radius, assuming a binary separation of 
0.1 AU. Given that not much is known about the binary 
fraction in the Galactic center, this is fully consistent 
with the observed number of 15 S stars. 

5. DISCUSSION 

The mechanism described above naturally explains the 
reported mysteries of the S stars: The disruption of stel- 
lar binaries as tight as 0.1-1 AU can leave a single star 
orbiting the central SMBH a factor of 10 closer than the 
original binary. This is just the distance ratio of S stars 
and stellar disks observed in the center of our Galaxy. 

The orientation of the disks rapidly changes due to 
precession. Furthermore, the change of orientation of 
a binary's orbit is a natural consequence of the Kozai 
effect. Once a binary is disrupted and leaves a single star 
on a much closer orbit, this orbit is decoupled from the 
originating stellar disk's motion, as both the precession 
frequency and the Kozai oscillation period are functions 
of central distance, explaining the widely varying orbital 
orientations observed. 

The eccentricities created for the tightly bound stars 
after tida l disruption are as high a s 0.99. Resonant re- 
laxation (jRauch &: Tremaind |1996[ ) and Kozai interac- 
tions between the S stars are expected to turn these ex- 
tremely high eccentricities into an isotropic distribution 
on a short timescale, just as is observed for the S stars. 

While all of the S stars observed so far are main- 
sequence B-type stars, observations of the stellar disks 
in the Galactic center mostly revealed OB giants and 
Wolf-Rayet stars. But even for the unusually flat mass 



function suggested (jPaumard et al.|[2006[ ). one would ex- 
pect a factor of 5 more B-type main-sequence stars in the 
disks than O stars. Most of these stars have not been 
observed yet, mainly due to crowding of the observed re- 
gion, limited spatial reso lution, and strong nebular emis- 
sion (Paumard et al.ll2006) . Our computation assuming 
equal- mass binaries predicts only one S star above 20 Mq 
per 18 B-type S stars. The lack of observed 0-type S 
stars is therefore not statistically significant. The num- 
ber of expected massive S stars would drop further for 
random pairing of binary component masses, since dis- 
ruption of a high mass ratio binary will only moderately 
change the orbit of the massive member. Furthermore, 
the massive O stars might not form in as tight binaries 
as are required to produce S stars. Hence, it is legiti- 
mate to assume that both S stars and disk stars formed 
in the same environment, thus providing a solution to 
the "paradox of youth" . 

In this Letter, we studied the evolution of two iso- 
lated stellar disks orbiting a SMBH. In a forthcoming 
publication, we will analyze the effects of other con- 
tributions of a perturbing potent ial, such as a stel- 
lar cusp fe.g.. iSchodel et aLri2007D. a massive circum- 
nuclear disk (IChristoph er et al]|2005f) . an intermediate- 
mass black hole and/or a star cluster (probably IRS 13; 
iMaillard et al. 2004). Preliminary results indicate, e.g., 
that the presence of a cusp of stellar black holes (e.g., 
iFreitag et al.l 120061 ) does not prevent the dynamical for- 
mation of S stars as described in this Letter. 

We thank Ladislav Subr for helpful discussions. This 
work was supported by the German Research Foundation 
(DFG) through the priority program 1177 "Witnesses 
of Cosmic History: Formation and Evolution of Black 
Holes, Galaxies and Their Environment." 
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